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Abstract 

In this paper we consider two different geometry of quasi one-dimensional semi- 
conductors and calculate their exchange-correlation induced bandgap renormaliza- 
tion (BGR) as a function of the electron-hole plasma density and quantum wire 
width. Based on different fabrication scheme, we define suitable external con- 
finement potential and then leading-order GW dynamical screening approximation 
is used in the calculation by treating electron-electron Coulomb interaction and 
electron-optical phonon interaction. Using a numerical scheme, screened Coulomb 
potential, probability of different states, profile of charge density and the values of 
the renormalized gap energy are calculated and the effects of variation of confine- 
ment potential width and temperature are studied. 
PACS: 73.20.Dx; 71.35.Ee; 71.45. Cm 
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1 Introduction 



A highly dense electron-hole plasma can be generated in a wide variety of semiconduc- 
tors by optical pumping. The band structure and the optical properties of highly excited 
semiconductors differ from those calculated for non-interacting electron-hole pairs due 
to many-body exchange-correlation effects arising from the electron-hole plasma [1-3]. 
One of the important many-body effect in high density electron-hole plasma is a density- 
dependent rcnormalization of the fundamental band gap of the semiconductor, which 
causes an increasing absorption in the spectral region below the lowest exciton resonance. 
The exchange-correlation correction of the fundamental band gap due to the presence of 
free carriers (electrons in the conduction band and holes in the valence band) in the sys- 
tem is referred to as the band gap rcnormalization (BGR) effect. Optical nonlinearities, 
which are strongly influenced by Coulomb interaction in the electron-hole plasma, are 
typically associated with the bandgap rcnormalization phenomenon. 
The band gap re-normalization has been widely studied in bulk and quasi-two dimensional 
(quantum well) semiconductors [4-6]. In recent years, quasi-one dimensional semiconduc- 
tor quantum wires (QW) have been fabricated in variety of geometric shape with atomic 
scale deflnition, and QW optical properties have been studied for their potential device 
applications such as semiconductor lasers [7- 9]. There has, however, been httle work on 
the BGR in QW systems, both experimentally and theoretically. Recently different ge- 
ometries of quantum wires, such as rectangular, V-shaped and T-shaped quantum wires 
have been fabricated and studied, and various experimental technics for fabrication and 
growth of these structures have been developed [2, 3]. Square quantum well wires have 
been studied by Hu and Das Sarma[10]. They have calculated the value of the band gap 
re-normalization for this case in GW approximation. Band gap renormalization in pho- 
toexcited semiconductor quantum wire in GW approximation has been studied by Hwang 
and Das Sarma[ll]. Rinaldi and Cingolani have studied optical properties of ID quantum 
structures specially the case of V-shaped quantum wire[12]. They have considered con- 
finement potential of the form V{y) = — ^^^^f^^^y^ but they have not calculated the band 
gap re-normalization with this confinement potential theoretically. Bener and Haug have 
considered plasma-density dependence of the optical spectra for quasi-one-dimensional 
quantum well wires[13]. Tanatar has studied band gap re-normalization in quasi-one di- 
mensional systems in simple plasmon-polc (quasi-static) approximation[14]. On the other 
hand, T-shaped quantum wires recently have been considered by some authors. For ex- 
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ample Sedlmaier and his coworkers have studied band gap re-normahzation of modulation 
doped T-shaped quantum wires. They have presented self-consistent electronic structure 
calculations for this device. These calculations show a band gap re-normalization which, 
when corrected for excitonic energy and its screening, are largely insensitive to the ID 
electronic density [15]. Stopa, using density functional theory, has calculated the elec- 
tronic structure of a modulation doped and gated T-shaped quantum wire. He also has 
calculated the band gap re-normalization as a function of the density of conduction band 
electrons[16]. Lin, Chen and Chuu have found the dependence of the bound states of 
L-shaped and T-shaped quantum wires to some asymmetric parameter in an inhomo- 
geneous magnetic fields[17]. Sedlmaier and his coworkers have calculated the band-gap 
re-normalization of modulation doped quantum wires by considering the photolumines- 
cence spectra as a function of the one dimensional density specially in T-shaped quantum 
wires[18]. Nozari and Madadi, recently have proposed a theoretical framework for calcu- 
lation of band gap re-normalization in V-shaped quantum wires. They have calculated 
numerically the value of this re-normalization in various temperatures and carrier density 
and have found exact solutions for some geometry of this type of quantum wires in a 
simplified random phase approximation[19]. So far most of calculations have been done 
in the static screening approximation or in the simple plasmon-pole approximation, which 
is a simplified version of the random-phase approximation (RPA). The plasmon-pole ap- 
proximation consists of ignoring the weight in the single particle excitations and assuming 
that all free carrier contributions to the dynamical dielectric function lies at the effective 
plasma frequency. The advantages of the plasmon-pole approximation are its mathemat- 
ical simplicity and simple physical meaning. However, a certain degree of arbitrariness 
in the choice of the effective plasmon pole parameters leads to considerable difficulties 
in applying the theory to semiconductors with complex band structures. In this paper, 
we calculate the BGR of the V-shaped and T-shaped quantum wire structures based on 
the RPA dynamical screening (GW) scheme by taking into account special mathematical 
definitions for confinement potentials which differ considerably from existing literature 
specially [12] and [19]. In doing so, the full frequency dependent dielectric response in the 
two component one-dimensional electron-hole plasma is considered. 
The structure of the paper is as follows: In section 2 we provide formal theory of bandgap 
renormalization in a general quasi one-dimensional structure using dynamical random 
phase approximation. Section 3 considers V-shaped confinement potential. First a suit- 
able mathematical definition for this form of confinement potential is given. Then a 
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numerical scheme for solutions of exact equations is provided and various quantities are 
calculated numerically. Section 4 considers the same analysis for T-shaped confinement 
potential. Conclusions and discussion are given in section 5. 

2 BGR in Dynamical Random Phase Approximation 

The exchange-correlation induced correction of the fundamental band gap due to the 
presence of free carriers (electrons in the conduction band and holes in the valence band) 
in the system is referred to as the band gap re-normahzation effect. The BGR is given by 
the sum of the self-energies for electrons and holes at band edges 

A = i?e(Se(0, 0)) + Re{Eh{0, 0)), (1) 

where Re stands for real part. The total electronic self-energy within the leading order 
effective dynamical interaction (GW) in a two-component electron-hole plasma is 

Ee(^, c.) = z y - y —Go{k -q^u-u ) (2) 

where Go is the Green's function for noninteracting electron gas, Vg = Vc + Vph the total 
effective interaction and e{q,uj) = 1 — Vs{q,uj)Ilo{q,uj) the effective dynamical dielectric 
function. V^,, V^/^ and IlQ{q,uj) arc direct Coulomb interaction, the longitudinal optical 
phonon mediated electron-electron interaction, and irreducible polarizability, respectively. 
The self-energy can be separated into the frequency-independent exchange and correlation 
parts 

Si(A;,a;) = Sf"(A;) + Sf°'-(A;,a;), i = e,h (3) 

where, 

^f^{k)^-J^^^n^{k + qmq). (4) 

Here, ni?(/c + q) = 9{kF — \k + q\) is the Fermi function at T = and 9 is Heaviside step 
function. In GW approximation, the Tf(°'^{k,oj) can be written in the following form, 

Ef-(A;,a;) = Ef-(^,a;) + T.f\Ku), (5) 

where 



dq duj' 



oo 27r J-oo 27r \{^k+q — oj) — ioj' 



e{q,iuj') 



- 1 



(6) 
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and 



where $,{k) is defined as 



0{UJ - ^k+q) - ^{-^k+q) 



Vs{q) 



- 1 



27.2 



2m* 



(7) 



(8) 



and m* is effective electron mass[ll]. 

In general the effective quasi-one dimensional Coulomb interaction in the framework of 
dynamical random phase approximation is given by 



Vs{k) 
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+00 

-00 ^^[[x - x'Y + {y - y'Y]^ 
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+00 
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+00 



dy'i.{k,y-y')\ip{y)\'\ip{y')\\ 



(9) 



In this equation ip{y) is transverse wave function and ^{k^y — y') = ^Ko[\k{y — y')|] 
where Ko{x) is zeroth order modified Bessel function of second kind. As we will show, in 
our situation the problem is effectively two dimensional (quantum well wire structure). 



3 BGR in V-Shaped Quantum Wires 

Using formalism of the last section, here we calculate density of states(profile of charges 
distribution) and band gap re- normalization in the geometry of V-grooved quantum quan- 
tum wire. There are many possibilities to define mathematically the geometry of V- 
grooved confinement potential[3,12]. These possibilities have originated from different 
fabrication and growing scheme. A possible and in some sense, technically suitable form 
of V-shaped confinement potential can be written as. 



V{x,y) 



{—Wx < X < Wx 
ax < y < ax + Wy 
—ax < y < —ax + 
00 elsewhere 

V 



(10) 
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For simplicity in numerical calculations, we define re-scaled quantities x — y — ^ 
and a — for x, y and a respectively. So the re-scaled confinement potential becomes, 



V{x,y) 



-l<x <1 
if < ax < y < ax + 1 

—ax < y < —ax + 1 
oo elsewhere. 



With this form of re-scaled confinement potential, one can compute band gap narrowing in 
V-grooved quantum wire. To do this end, one should calculate effective two-dimensional 
transverse wave function. According to above geometry, z-direction is without confine- 
ment, so carriers move in this direction freely. Now the Schrodinger equation for carrier 
in two direction of confinement is. 



This can be written, using re-scaled quantities, as 

1 d"^ 1 9^ 

(-7^)( — ^ + — ^)'^(^'^) + V{x,y)(i){x,y) = E(p{x,y). 



(12) 



2m dx'^ Wy r>^''^ ' 



Defining, E — V — ^ and c 



(13) 



'-, we find finally 



dx"^ dy^ 



(14) 



This equation with potential as (10), can be solved analytically. The screened coulomb 
potential now is given by, 

2e2 



Vs{k) = ^J dxdy J dx'dy'Ko\kArMx,y)\^\^{x',y')\\ 



(15) 



or 



2e2 



Vs{k) = \=^wlwl J dx'dy'J di'dy'Ko\kArMx,y)\'\ip{x',y')\', (16) 



where A is re-scaling factor equal to 10~^^/m^ and kAr = ky {x — Xq)"^ + (y — yoY- Since 
we want to solve the integral of equation (16), and this integral can not be solved analyti- 
cally, wc need to compute it numerically and therefore the ground state energy and ground 
state wave function of Schrodinger equation should be calculated numerically. Therefore, 
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we first provide a suitable numerical scheme for calculations. 

Basically, we have to solve 2D-Schrodinger equation (13) for potential as given in equa- 
tion (10), which is not, theoretically, solvable yet. Therefore, we are going to calculate 
the ground state wave function numerically. We use the usual finite difference algorithm 
to solve this eigenvalue problem[20]. First, we do define the re-scaled values for x,y,E 
and V respectively as ^"^J"'^ and ^!!^^. Then, we do discritize the x,y axes's to 

discritized space dx — 0.1 and dy — 0.1, therefore the equation (13) will now be changed 
to form of the matrix equation as: 

Hep = E(f), (17) 

where 7i, and are the Hamiltonian matrix and state wave function array respectively 
and are defined as: 



= 



/ 



(p{xi+i,yj) 



\ 



J 



<P[Xi,yj+i) 
(j){xi+i,yj+i) 



(18) 
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and 
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2(^ + i^) + n^i,%-) 
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— c 

dx2 



^ dx^ dy 



+ V{xi+i,yj) 








V 







1 

dp 











(19) 

The error of computation of (f)[x, y) is in the order of 0{dx^). We diagonahze the matrix 
and calculate the ground state wave function numerically. The screened potential is 
calculated by ground state wave function as a function of k. Figure 1 shows the calculated 
screened potential versus the wave vector for different a. In this figure, Vs{k) is normalized 
by ^ and the k is normalized to kwy. Now, we are able to calculate the self energies of 
electrons and holes. For this goal, first wc define the rc-scalcd f3e and respectively as 
2^^2 = f^^w^ ^"^^ f^eP^ where Wy is width of quantum well in y direction in Nanometer. 
For holes, we also define the re-scaled Ph and as (3h = Pe{^) and = Aie(^) 

h h 

respectively. In all computations in this paper, we have assumed that the ratio ^ is 
equal to 0.3 and m* ~ O.OGTmg. By calculating the self energies of holes and electrons, 
wc compute the band gap rc-normalized energy for V-groovcd confinement potential in 
dynamical random phase(GW) approximation. We consider the case of T = 0. The 
result is shown in figure 2 which shows the values of bandgap renormalization in terms 
of the width of the V-shaped confinement potential. In all figures, we have assumed that 
chemical potentials for electron is ImeV and Eg is normalized by 2e^/eo. 
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4 BGR in T-Shaped Quantum Wires 



Now we want to do the same calculations for geometry of T-shaped quantum wires. It is 
evident that one can consider the external T-shaped confinement potential as follows, 



V{x,y) = < 






oo 



— oo < X < oo. 



'-^<y<oo, 



--f^<x<^ 



(20) 



elsewhere. 



Note that we want to simulate actual situation with mathematical functions. The actual 
problem is a quantum wire which is a quasi-one dimensional structure but its external 
confinement potential should be effectively two dimensional to provide two direction of 
confinement. Therefore, this form of confinement potential is considered effectively two 
dimensional. To proceed the same algorithm as preceding section, we first re-scale some 
quantities. First of all the confinement potential can be written as. 



V ^ i 






oo 



-oo < X < oo. 



^ < y < oo. 



elsewhere. 



(21) 



The re-scaled Schrodinger equation becomes 



-( 



dx'^ 



+ a 



)(f){x, y) + — -2-^ — 0(x, y) 



dy^ 



h'wl 



2mE 



(22) 



where now a is ratio of well width potential in x direction to the y direction, i.e. a — 
The same numerical scheme has been used to calculate various quantities. By diagonaliza- 
tion of the Hamiltonian matrix, we calculate the ground state wave function numerically. 
Using ground state wave function, one can calculate the concentration of carriers in quan- 
tum wire structure and therefore density of probability for these carriers. This procedure 
is as previous section. One can show that by decreasing a, the wave function of carriers 
will be localized in x directions while by increasing a the wave function will be more 
localized in y directions of the well relative to x direction. From physical grounds, this is 
reasonable to say that probability density and charge distribution are function of external 
confinement potential. 

Now we should calculate the screened Coulomb potential. In our situation the screened 
Coulomb potential should be considered more carefully, since now the confinement poten- 
tial of quantum well wire is effectively two dimensional. It is important to note that 
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this means that quantum mechanical problem under consideration is effectively two- 
dimensional. For these reasons, the screened Coulomb potential can be written as 

dxdy / dx'dy'u{k,p- p')\ip{x,y)\^\ip{x',y')f (23) 

where p = \xi + yj\ is the radius polar coordinate in plane and Lp{x,y) is transverse 
wave function and v{k^p — p') = ^KQ[\k{p — p')\] where Kq(x) is zeroth order modified 
Bessel function of second kind. With this point in mind the screened Coulomb potential 
(16) is calculated by two dimensional ground state wave function. We have calculated 
the screened Coulomb potential, as a function of k in figure (3). In this figure, Vs{k) is 
normalized by ^ and the k is normalized to 1/Wy and the figure shows that screened 
potential is going to zero at lower wave length, and it is infinite at higher wave length 
for different a ratio. The screened potential for a constant wave vector of carriers will 
increase by increasing a. This indicates that screening effects for carriers are affected by 
the ratio of width of well potentials. 

Now to calculate the self energies of electrons and holes first we define the re-scaled Pe 
and /Xg respectively as before. By calculating the self energies of holes and electrons, we 
compute the band gap energy renormalization for different temperatures and width of the 
quantum wells for T-shaped potential. Figure (4) shows the relative difference of band 
gap energy renormalization for T-shaped potential versus the width of quantum well in y 
directions. Increasing temperature leads to more relative renormalization of gap energy 
and this is in agreement with the experimental results [1,1 3]. Also it is evident that in our 
situation the relative width of the well, i.e. a has significant effect on the value of band 
gap renormalization. Also our analysis shows a temperature dependence of gap energy 
renormalization. Generally, under high optical excitation, the band gap for 2D and bulk 
systems is found to decrease with plasma density due to exchange-correlation effects. The 
observed band gaps are typically normalized by ~ 20meV within the range of plasma 
densities of interest which arise from the conduction band electrons and valence band 
holes. Since Coulomb screening decreases with reduction of semiconductor dimensions, 
typical values of band gap renormalization in the case of bulk semiconductor is more than 
gap renormalization in quantum well structures. In the same manner typical band gap 
renormalization for quantum wires is less than that of quantum wells. Our numerical 
results, when are compared with experimental data[13], show the good agreement with 
experiments at least in general behavior. 



10 



5 Results and Discussions 



In this paper, a numerical approach has been proposed to calculate band gap renormaliza- 
tion in a V-shaped and T-shaped quantum wires. Quantum field theoretical random phase 
approximation within the leading-order screening has been used to calculate screened 
Coulomb potential in terms of carrier density and width of the confinement potentials. 
The essential question is about the form of confinement potentials. As we have argued 
the quantum mechanical problem under consideration is effectively two dimensional. In 
other words since we have two dimension of confinement, simulation of external confine- 
ment with mathematical functions, leads to an effectively two dimensional confinement 
potential which requires solution of two dimensional Schrodinger equation. Based on this 
viewpoint we have shown that using two dimensional Cartesian and polar coordinates for 
V-shaped and T-shaped quantum wire respectively, one can use the formalism of Hu and 
Das Sarma with re-definition of some quantities. Our numerical calculations show that 
the distribution of probability for carriers concentration is a function of the ratio of the 
well width in two dimensions. This is natural since probability distributions is a function 
of confinement potential. In other words, the geometric shape of external confinement 
itself restricts the shape of wave functions. As the figures show, by decreasing a, the wave 
function of carriers will be localized in x directions while by increasing a the wave function 
will be more localized in y directions of the well relative to x direction. It is important to 
note that by re-scaling procedure which we have considered, we have fixed the geometry 
of the wires but now Schrodinger equation becomes re-scaled via the presence of ratio 
a. Actually one should consider the possibility for changing the geometry also. But our 
investigation show that the main physical results do not change considerably. Calculation 
of screened Coulomb interaction shows that the screened potential is going to zero at lower 
wave length, and it is infinite at higher wave length for different a ratio. The screened po- 
tential for a constant wave vector of carriers will increase by increasing a. This indicates 
that screening effects for carriers are affected by the ratio of width of well potentials. Also 
calculation of band gap renormalization shows that increasing temperature leads to more 
relative renormalization of gap energy and this is in agreement with the experimental 
results. Also it is evident that in our situation the relative width of the well, i.e. a has 
significant effect on the value of band gap renormalization. This is not surprising since the 
screening effect itself is dependent on the shape of confinement potential. Since Coulomb 
screening decreases with reduction of semiconductor dimensions, typical values of band 
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gap renormalization in the case of bulk semiconductor is more than quantum well struc- 
tures. In the same manner typical band gap renormalization for quantum wires is less 
than that of quantum wells. The relative values of band gap renormalization in various 
fabrication pattern(i.e. V-shaped, T-shaped and rectangular) quantum wire structure is 
dependent to the shape and depth of the potential well. Our numerical results, when 
are compared with experimental data[13], show the good agreement with experiments at 
least in general behavior. 
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Figure 2: Renormalization of gap energy calculated numerically for V-shaped potential versus 
the width of quantum well wire Wy. In all figures, we have assumed that chemical potentials for 
electron is ImeV. The AEg is normalized by 2e^/eo. 
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Figure 3: The calculated screened Coulomb potential versus the wave vector. The screened 
potential is normalized by ^ and the k is normalized to 1/wy. 
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Figure 4: This figure shows the difference of band gap renormalized energy for T-shaped 
confinement potential versus the width of quantum well in y directions for different temperatures 
and different ratio of the well width in two dimensions, a. In all figures, we have assumed that 
the chemical potentials for electron is ImeV. The AE'^ is normalized by 2e^/eQ. 
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